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Control  Allocation  for  Overactuated  Systems 

Michael  W.  Oppenheimer  and  David  B.  Doman 


I.  Abstract 

Much  emphasis  has  been  placed  on  overactuated  systems 
for  air  vehicles.  Overactuating  an  air  vehicle  provides  a 
certain  amount  of  redundancy  for  the  flight  control  system, 
thus  potentially  allowing  for  recovery  from  off-nominal 
conditions.  Due  to  this  redundancy,  control  allocation  algo¬ 
rithms  are  typically  utilized  to  compute  a  unique  solution 
to  the  overactuated  problem.  Control  allocators  compute 
the  commands  that  are  applied  to  the  actuators  so  that  a 
certain  set  of  forces  or  moments  are  generated  by  the  control 
effectors.  Usually,  control  allocation  problems  are  formulated 
as  optimization  problems  so  that  all  of  the  available  degrees 
of  freedom  can  be  utilized  and,  when  sufficient  control  power 
exists,  secondary  objectives  can  be  achieved.  In  this  work,  a 
survey  of  control  allocation  techniques  will  be  given. 

II.  Introduction 

Conventional  aircraft  utilize  an  elevator  for  pitch  control, 
ailerons  for  roll  control,  and  a  rudder  for  yaw  control. 
As  aircraft  design  advances,  more  control  effectors  (some 
unconventional)  are  being  placed  on  the  vehicles.  In  some 
cases,  certain  control  effectors  may  be  able  to  exert  sig¬ 
nificant  influence  upon  multiple  axes.  When  a  system  is 
equipped  with  more  effectors  than  axes  to  control,  the  system 
may  be  overactuated  or  redundant.  The  allocation,  blending, 
or  mixing  of  these  control  effectors  to  achieve  some  desired 
objectives  is  the  control  allocation  problem. 

Due  to  overactuation  and  coupling  of  control  surface 
effects,  it  is  difficult  to  determine  an  appropriate  method 
of  how  to  translate  a  flight  control  command  into  a  control 
surface  command.  In  addition,  rate  and  position  limits  of  the 
control  surfaces  must  be  considered  in  order  to  achieve  a 
realistic  solution.  Not  only  is  the  mixing  of  control  surface 
effects  critical,  it  is  also  desired  to  enable  the  aircraft  to 
recover  from  off-nominal  conditions,  such  as  a  failed  control 
surface,  when  physically  possible.  In  reconfigurable  control 
systems,  a  control  allocation  algorithm  is  needed  to  perform 
automatic  distribution  of  the  control  power  requests  among 
a  large  number  of  control  effectors,  while  still  obeying  the 
rate  and  position  limits  of  the  actuators  and  to  potentially 
allow  recovery  from  off-nominal  conditions. 

Some  of  the  simplest  control  allocation  techniques  are 
explicit  ganging,  pseudo  control,  pseudo  inverse,  and  daisy 
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chaining.  Unfortunately,  each  suffers  from  difficulty  in  guar¬ 
anteeing  that  rate  and  position  limits  will  not  be  violated 
and  some  can  be  difficult  to  apply  due  to  the  need  to 
derive  a  control  mixing  law  a  priori.  Another  control  allo¬ 
cation  method,  called  direct  allocation  [1],  finds  the  con¬ 
trol  vector  that  results  in  the  best  approximation  of  the 
command  vector  in  a  given  direction.  Unconstrained  least 
squares  control  allocation  methods,  that  account  for  rate 
and  position  limits,  through  the  use  of  penalty  functions, 
have  also  been  developed  [2].  One  of  the  first  instances 
of  linear  programming  based  control  allocators  was  from 
Paradiso  [3].  In  this  work,  Paradiso  developed  a  selection 
procedure  for  determining  actuator  positions  that  was  based 
on  linear  programming  and  limited  actuator  authority.  More 
recently,  the  control  allocation  paradigm  has  been  posed  as 
a  constrained  optimization  problem  [4].  In  this  work,  the 
control  allocation  problem  was  split  into  two  sub-problems. 
The  first  was  the  error  minimization  part,  which  attempts 
to  find  the  control  vector,  such  that  the  control  effector 
induced  moments  or  accelerations  match  the  desired  mo¬ 
ments  or  accelerations.  If  multiple  solutions  exist  to  the 
error  minimization  problem,  the  second  problem  attempts 
to  find  a  unique  solution  by  driving  the  control  vector  to 
a  preference  vector  and  optimizing  a  secondary  objective. 
The  linear  control  allocation  problem  has  been  extended  to 
an  affine  problem  [5]  to  account  for  nonlinearities  in  the 
moment-deflection  curves.  Quadratic  programming  has  also 
been  used  in  the  past  [6].  An  excellent  paper  discussing 
control  allocation,  by  Bodson  [7],  provides  a  glimpse  into 
numerous  control  allocation  techniques. 

This  paper  will  cover  some  of  the  more  popular  linear 
control  allocation  techniques.  The  objective  is  to  present  the 
methods  in  sufficient  detail  that  the  readers  can  utilize  the 
techniques. 

III.  Linear  Control  Allocation 

The  linear  control  allocation  problem  can  be  posed  as 
follows:  find  the  control  vector,  <5  £  Rn,  such  that 

B<5  =  Ades  (1) 

subject  to 

&min  —  &  —  &max 

k  k  \^) 

"  _  "max 

where  B  £  Rrax"  is  a  control  effectiveness  matrix,  the 
lower  and  upper  position  limits  are  defined  by  dmin  £  R" 
and  6max  £  R",  respectively,  5  £  R"  are  the  control 
rates,  Smax  £  R"  are  the  maximum  control  rates,  Ades  are 
the  desired  moments  or  accelerations  (typically  for  inner- 
loop  control  laws,  d,/,.s  £  R3),  n  is  the  number  of  control 
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effectors,  and  m  is  the  number  of  axes  to  control.  Equation  2 
provides  the  position  and  rate  limits  for  the  control  effectors. 
In  a  digital  computer  implementation,  the  rate  limits  are 
converted  to  effective  position  limits.  The  combined  limits 
become  the  most  restrictive  of  the  rate  or  position  limits  and 
are  specified  as 

8  <  S  <  8  (3) 


where 


8  = 
8  = 


mm  dmax 


8  +  A  tSn 


max  (  5 min .  8  —  A tdr 


(4) 


Here,  8  £  R",  <5  £  R",  and  8,  8  are  the  most  restrictive 
upper  and  lower  control  effector  limits,  respectively. 

A  necessary  condition  for  a  system  to  be  overactuated  is 
the  number  of  columns  of  B,  n,  must  be  greater  than  the 
number  of  rows  of  B,  m.  The  true  test  of  overactuation  is  that 
the  number  of  linearly  independent  columns  of  B  be  greater 
than  the  number  of  rows  of  B.  For  inner-loop  control  laws, 
the  B,  or  control  effectiveness,  matrix  typically  becomes 


B  = 


dL_  dL_  . 
<9<5i  <9<52 
3M  3M 
<9<5i  <9<52 
8N_dN_  . 
<9<5i  <9<52 


dL 

dSn 

dM 

d5n 

dN 

dSn 


(5) 


where  L,  M,  and  N  are  the  rolling,  pitching,  and  yawing 
moments,  respectively. 

The  control  allocation  problem  can  be  illustrated  with  a 
simple  example.  Consider  the  following  problem:  find  <5i  and 
62  such  that 

ddes  =  3<5i  +  62  (6) 


In  this  case,  there  is  one  objective  function,  ddes  €E  R1  and 
two  controls.  Using  the  form  shown  in  Eq.  1  yields 

B  =  [  3  1  ]  S  =  [  <fi  62  ] T  (7) 

Let  the  most  restrictive  limits  be 


-1  <  <5i  <  1  -1  <  62  <  1  (8) 


Let  ddes  =  2.  Then,  this  problem  could  be  solved  in  a 
number  of  ways,  for  example,  Si  =  0  and  6 2  =  2,  or 
(5i  =  |  and  82  =  0,  or  <5i  =  —3  and  82  =  11,  or 
many  others.  Graphically,  this  situation  becomes  as  shown 
in  Figure  1 .  The  control  allocation  problem  is  to  find  <5i ,  82 
such  that  ddes  =  3<5i  -EA2  and  the  constraints  are  not  violated. 
The  solution  is  the  intersection  of  the  hyperspaces  of  the 
constraints  and  the  equation  ddes  =  3<5i  +  82-  For  ddes  =  2 
or  3,  multiple  solutions  exist.  For  ddes  —  4,  only  one  solution 
exists,  while  for  ddes  =  5,  no  solutions  exist.  When  only  one 
solution  exists,  we  simply  select  that  solution.  When  multiple 
solutions  exist,  a  method  to  pick  one  is  necessary.  When  no 
solutions  exist,  a  method  to  minimize  the  error,  between  ddes 
and  35i+(>2  is  required,  in  some  logical  sense.  This  illustrates 
the  control  allocation  problem. 

Equations  1,  3,  and  4  define  the  linear  control  allocation 
problem.  The  objective  now  is  to  determine  methods  which 


Fig.  1.  Example  of  Control  Allocation. 


allow  computation  of  the  control  effector  vector  <5,  while  pos¬ 
sibly  taking  into  account  effector  rate  and  position  limits.  The 
following  discussion  looks  at  different  methods,  which  can 
be  utilized  to  either  reduce  the  dimension  of  the  overactuated 
system  to  the  point  that  a  square  allocation  problem  results 
or  to  directly  solve  the  overactuated  linear  control  allocation 
problem. 

A.  Explicit  Ganging 

In  this  approach,  an  a  priori  method  is  used  to  combine  or 
gang  the  effectors  to  produce  a  single  effective  control  from 
multiple  devices.  Note  that  historically,  ganging  was  done 
with  cables,  pulleys,  or  other  mechanical  means.  On  modern 
aircraft,  fly-by-wire  is  used  and  the  ganging  is  performed 
in  software.  The  goal  is  to  find  a  matrix  G,  that  relates  the 
pseudo  controls,  Spseudo ,  to  the  actual  controls,  S  such  that 

S  =  G  SpSeudo  (9) 

The  vector  Spseudo  are  the  pseudo  controls,  so  named  be¬ 
cause  some  or  all  components  of  this  vector  are  not  physical 
control  surfaces  on  the  vehicle.  For  example,  consider  a 
vehicle  that  has  left  and  right  ailerons  for  roll  control, 
SaL  £  R1  and  8aR  £  R1,  a  single  rudder  for  yaw  control, 
Sr  £  R1,  and  a  single  elevator  for  pitch  control,  8e  £  R1.  A 
priori,  a  ganging  law  is  constructed  to  produce  a  single  roll 
control  device.  One  possibility  is  to  let 

<5a  =  0.5  (SaL  +  SaR)  (10) 

where  5a  £  R1  is  the  single  roll  control  device.  Therefore, 
the  full  ganging  law  becomes 


$aR 

'100' 
10  0 

'  <5a  ' 

£ 

5e 

0  10 

8r 

0  0  1 

or 

3  —  G$  pseudo  (12) 

Here,  the  reason  for  the  term  pseudo  controls  becomes  clear 
because  one  of  the  elements  of  Spseud0,  namely,  5a,  is  not 
a  physical  control,  instead,  it  is  a  linear  combination  of  two 
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physical  control  effectors,  SaL  and  S„  R .  Then,  the  control 
allocation  problem  becomes  find  Spseud0 ,  such  that 

B<5  =  ddes  RGSpSeudo  —  drfes  (13) 

Solving  this  allocation  problem  for  Spseud0  and  using  Eq.  12 
yields  the  physical  control  effector  commands.  Typically, 
explicit  ganging  is  used  when  it  is  obvious  how  to  combine 
the  redundant  control  effectors. 

It  is  important  to  point  out  that  this  method  can  be  used  to 
reduce  the  control  space  dimension  of  an  overactuated  sys¬ 
tem.  As  previously  mentioned,  the  inner-loop  aircraft  control 
law  typically  contains  three  objective  functions,  namely  that 
the  moments  produced  by  the  controls  are  equal  to  a  set  of 
desired  moments  (d,i,s)-  In  the  example  problem,  d  £  R4 , 
so  that  B<5  =  &des  is  a  nonsquare  allocation  problem.  After 
employing  the  explicit  ganging  methodology,  Spseudo  £  R3 
so  that  BG Spseudo  =  d des  is  a  square  allocation  problem. 
Hence,  the  dimension  of  the  control  space  is  reduced  from 
4  to  3. 

B.  Pseudo  Inverse 

The  pseudo  inverse  method  is  a  constrained  optimization 
technique  that  requires  a  pseudo  inversion  of  the  generally 
nonsquare  B  matrix.  The  pseudo  inverse  solution  is  the  two- 
norm  solution  to  the  control  allocation  problem  and  can  be 
formulated  as  follows: 

min  J  =  min  -(6  +  c)TW(5  +  c)  (14) 

s  s  2 

subject  to 

B  S  =  ddes  (15) 


Substituting  Equation  20  into  Equation  17  produces 

W£  =  -Wc  +  BT(BW~1BT)-1[ddes  +Bc]  (21) 
Simplifying  Equation  21  gives  the  desired  result 

^  =  — c  +  W_1BT(BW_1BT)_1  [ddes  +  Be]  ^ 
= -c  +  B»[ddes+Bc] 

where  B'  is  the  pseudo  inverse  of  B  and  the  superscript  (**) 
indicates  a  pseudo  inverse  operation.  Equation  22  gives  the 
pseudo  inverse  solution.  It  should  be  noted  that  if  an  effector 
is  offset,  two  items  must  be  taken  into  account,  position  offset 
(— c)  and  the  moments  generated  by  the  offset  (Be).  For  the 
position  offset,  the  negative  of  the  locked  position  is  placed  in 
the  corresponding  entry  of  the  c  vector.  For  example,  assume 
there  are  four  controls  and  control  number  3  is  stuck  at  +5 
degrees.  Then,  the  c  vector  would  become 

c  =  [  0  0  -5  0  ]T  (23) 

The  weighting  matrix,  W,  can  be  selected  to  incorporate 
the  position  limits  of  the  control  effectors.  For  example,  a 
diagonal  element  of  W  can  be  selected  to  be  a  function  of 
the  corresponding  component  of  5,  so  that  the  weighting 
function  approaches  oo  as  the  control  approaches  a  physical 
limit.  There  are  no  guarantees  that  commands  to  the  control 
effectors  will  not  exceed  the  position  limits;  however,  in 
practice  the  method  is  effective  in  constraining  the  positions 
of  the  controls.  As  a  final  note,  this  method  can  be  useful  in 
generating  preference  vectors  for  more  complex  optimization 
based  methods  for  the  purpose  of  robustness  analysis. 

C.  Redistributed  Pseudo  Inverse 


where  W  £  Rraxn  is  a  weighting  matrix  and  c  €  R"  is 
an  offset  vector  used  to  represent  an  off-nominal  condition 
with  one  or  more  control  effectors.  To  solve  this  problem, 
first  find  the  Hamiltonian  (iT)  such  that 


H=-  (sTWS 

2  V 

+  £(B5  —  ddes) 


cTWS  +  STWc  +  cTWc 


(16) 


where  £  £  Rn  is  an  as  yet  undetermined  Lagrange  multiplier. 
Taking  the  partial  derivatives  of  H  with  respect  to  S  and  £, 
setting  these  expressions  equal  to  zero,  and  rearranging,  gives 


W  =  W^^(cTwf  +  iwc- 


(4B)T  =  0 


W<5  =  -  Wc  -  Bt£t 


and 


DH 

B 6  =  d des  =>  BW 


=  B(5  -  ddes  =  o 


(17) 


(18) 


‘W  S  =  d^s 

Substituting  Equation  17  into  Equation  18  yields 
BW-1[— Wc  —  Bt£t]  =  ddes 
Solving  for  £7  in  Equation  19  yields 

£T  =  -(BW_1BT)_1[ddes 


The  redistributed  pseudo  inverse  works  in  a  fashion  sim¬ 
ilar  to  the  pseudo  inverse.  The  difference  is  that  now  the 
process  is  iterative  and  position  saturated  control  effectors 
are  removed  from  subsequent  pseudo  inverse  solutions.  The 
first  step  is  to  solve  the  control  allocation  problem  using  the 
pseudo  inverse  solution  in  Eq.  22  with  c  initially  a  vector  of 
all  zeros.  If  no  controls  exceed  their  minimum  or  maximum 
position  limits,  then  the  process  stops  and  the  solution  from 
Eq.  22  is  used.  However,  if  one  control  saturates,  the  problem 
is  solved  again,  this  time  zeroing  out  the  column  of  the  B 
matrix  corresponding  to  the  saturated  control  and  placing 
the  negative  of  the  saturated  value  in  the  vector  c.  One 
must  be  careful  here.  When  saturation  occurs,  there  are  two 
B  matrices  in  Eq.  22,  one  for  the  pseudo  inverse  solution 


W  B  (BW  B)  and  one  for  the  offset  or  saturated 
contribution  Be.  When  zeroing  out  a  column  corresponding 
to  a  saturated  effector,  only  the  pseudo  inverse  B  matrix 
is  modified,  while  the  Be  term  uses  the  original  B  matrix. 
Consider  the  following  example  of  a  redistributed  pseudo 


Bel 


inverse  control  allocation  problem.  Let 

(19) 

B  = 

'  2 

1 

-2 

1 

-2  -1  ' 
-3  -2 

ddes  =[-51 

1  ]T 

2 

-2 

-1  -1 

(24) 

(20) 

-0.75  <  <5i,2,3,4  <  0.75 

(25) 
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In  this  example,  there  are  3  objective  functions  (3  rows  in 
B),  four  control  effectors  (4  columns  in  B),  each  control  has 
the  same  lower  and  upper  position  limits,  and  the  weighting 
matrix  is  W  =  I.  The  first  step  in  the  redistributed  pseudo 
inverse  solution  is  to  compute  Eq.  22.  The  results  are 

8  =  B^ddes  ==  [  0.55  0.23  0.5  —0.86  ]  ^  (26) 

where  the  c  vector  was  removed  since  it  currently  contributes 
nothing.  Since  r)'i  exceeds  it  limit,  force  hi  to  —0.75  and  zero 
out  the  fourth  column  of  the  B  matrix.  Now,  calculate  the 
new  pseudo  inverse  solution  using 

8  =  -c  +  B“ed  [ddes  -  Be]  =  [  .6875  .3125  .5  -.75  ]  T 

(27) 

where  Breai  is  the  reduced  B  matrix  formed  by  zeroing  out 
the  column  corresponding  to  the  saturated  control  effector. 
The  results  show  that  5 4  is  now  at  its  most  negative  limit,  as 
expected,  and  no  other  control  effectors  exceed  their  limit. 
The  iterations  stop  and  the  result  of  the  redistributed  pseudo 
inverse  calculation  is 

8  =[.6875  .3125  .5  —.75  ]  ^  (28) 

A  check  of  this  result  gives 

B  [  .6875  .3125  .5  -.75  ]r=[0.5  1  l]T  =  ddes 

(29) 

Thus,  the  calculated  control  settings  do  provide  the  desired 
commands.  In  this  example,  only  2  iterations  of  pseudo 
inverse  calculations  were  performed.  If  the  control  effector 
setting  computed  in  Eq.  27  had  one  or  more  controls  ex¬ 
ceeding  their  position  limits,  then  the  process  of  zeroing  out 
columns  in  the  B  matrix  would  continue  until  all  controls 
are  at  their  limits  (either  positive  or  negative)  or  a  feasible 
solution  is  found  (one  where  remaining  controls  do  not 
exceed  position  limits). 

D.  Daisy  Chaining 

The  daisy  chain  assumes  a  hierarchy  of  control  effectors. 
In  this  method,  when  one  control  or  a  group  of  controls 
saturates,  there  is  an  error  between  the  commanded  moments 
or  accelerations  and  those  produced  by  the  control  effectors. 
The  daisy  chain  method  would  then  utilize  another  control 
to  produce  the  required  moments  or  accelerations  that  are 
lacking  due  to  the  saturation  of  a  control  effector.  Figure  2 
shows  an  example  of  daisy  chain  allocation.  In  this  example, 
the  goal  is  to  produce  a  desired  pitch  acceleration,  given 
by  qdes-  There  are  three  controls  that  can  produce  pitching 
moment,  an  elevator  (Se),  a  bodyflap  {Sbf),  and  a  canard  ( Sc ). 
The  daisy  chain  procedure  works  as  follows:  the  primary 
control  effector,  5e  in  this  case,  is  commanded  to  produce  the 
desired  acceleration.  If  the  elevator  can  produce  this  accelera¬ 
tion,  then  nothing  else  happens  and  the  bodyflap  and  canard 
are  not  utilized.  However,  if  there  is  a  moment  deficiency 
between  the  acceleration  that  the  elevator  produces  and  the 
desired  acceleration,  the  control  effector,  which  is  second 
in  line,  namely  the  bodyflap,  is  commanded  to  produce  an 
acceleration  equivalent  to  the  acceleration  deficiency.  If  the 
bodyflap  can  produce  the  required  acceleration,  then  the  ca¬ 
nard  does  nothing.  However,  if  the  bodyflap  cannot  produce 


Fig.  2.  Example  of  Daisy  Chain  Allocation. 


the  required  acceleration,  then  the  canard  is  commanded  to 
produce  the  difference  between  the  commanded  acceleration 
and  the  accelerations  produced  by  the  elevator  and  bodyflap. 
This  method  is  valid  for  any  number  of  control  effectors. 

E.  Direct  Allocation 

The  direct  allocation  method,  by  Durham  [1],  is  a  con¬ 
strained  control  allocation  method  aimed  at  finding  a  real 
number,  p,  and  a  vector,  <5i,  such  that 

B61  =  pddes  (30) 

and 

8<8<8  (31) 

If  p  >  1,  let  8  =  ±8 1.  If  p  <  1,  let  (5  =  81.  In  order 
to  use  this  method,  an  attainable  moment  set  [9]  (AMS) 
must  be  established.  The  AMS  is  of  the  dimension  of  ddes 
and  for  linear  systems  of  the  form  of  Eq.  30,  consists  of 
a  convex  hull  with  planar  surfaces.  Physically,  p  represents 
how  much  a  control  power  demand  must  be  scaled  to  touch 
the  boundary  of  the  AMS.  When  p  <  1,  the  moment  demand 
lies  within  the  AMS  and  the  allocator  supplies  the  demand. 
When  p  >  1,  the  control  power  demand  exceeds  supply  and 
the  demand  is  scaled  back  to  touch  the  boundary  of  the  AMS. 
Algorithms  exist  for  generating  the  AMS  for  the  3  moment 
problem  [9]  and  the  direct  allocation  problem  extends  to 
nonlinear  problems  of  the  form  f{8)  =  pddes-  Unfortunately, 
the  construction  of  the  AMS  for  a  nonlinear  problem  can  be 
extremely  difficult. 

F.  Error  and  Control  Minimization 

The  objective  of  the  error  minimization,  control  law  com¬ 
mand  feasibility,  or  control  deficiency  problem  is  to  find  a 
vector  8 ,  given  a  matrix  B,  such  that 

min  J  =  ||B(S  —  ddes||  (32) 

s 

is  minimized,  subject  to 

8  <  8  <  8  (33) 

The  norm  used  would  depend  on  the  type  of  algorithm 
used  to  perform  the  minimization.  When  used  with  linear 
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programming  solvers,  the  error  minimization  problem  is 
specified  as 


minJ=  ||Btf  —  drfe»||i  (34) 
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subject  to  the  constraints  specified  in  Eq.  33.  This  can 
be  transformed  into  the  standard  linear  programming  prob¬ 
lem  [4] 


min  J  =  0 

s  L 


0  1  •  •  •  1 


S 

ss 


subject  to 


<5S 

0 

-s 

-6 

6 

> 

6 

— B<5  +  6S 

d des 

B(5  +  <5S 

&des 

(35) 


(36) 


where  S3  £  Rm  is  a  vector  of  slack  variables.  Note  that 
the  slack  variables  must  be  positive,  but  are  otherwise 
unconstrained.  Individually,  the  slack  variables  represent  how 
much  control  power  demand  exceeds  supply  in  any  given 
axis.  If  J  =  0,  then  the  control  law  command  is  feasible, 
otherwise  it  is  infeasible.  The  solution  to  Eqs.  35  and  36 
can  be  obtained  with  a  linear  programming  solver. 

The  control  minimization  or  control  sufficiency  problem 
is  a  secondary  optimization.  If  there  exists  sufficient  control 
authority  to  satisfy  Eq.  34,  that  is,  if  J  =  0,  then  a  secondary 
objective  may  be  achieved.  The  ability  to  do  this  is  a  direct 
result  of  the  overactuated  system,  in  that  multiple  solutions 
to  the  problem  may  exist  and  one  solution  may  be  preferred 
over  another.  The  control  minimization  is  posed  as  follows: 


min  J  =  Wlds 
8 

subject  to 


Ss 

0 

-6 

-s 

6 

> 

6 

—S  +  Ss 

1 

73 

s  +  s3 

1 

73 

1 _ 

B<5  =  ddes 


(37) 


(38) 


where  £  Rmxl,  Ss  £  Rm,  and  6P  £  R"  is  the  preferred 
control  effector  position  vector.  So,  the  first  requirement  is 
that  the  moment  or  acceleration  demand  is  met,  followed 
by  selecting  the  control  effector  positions  which  satisfy 
the  moment  demand  and  minimize  a  secondary  objective. 
Many  secondary  objectives  exist,  for  example,  minimum 
control  deflection,  minimum  wing  loading,  minimum  radar 
signature,  minimum  drag,  minimum  actuator  power,  and 
others.  A  few  of  these  secondary  objectives  are  discussed 
in  [4], 


G.  Minimum  2-Norm  for  Robustness  Analysis 

One  of  the  drawbacks  of  the  optimization  based  linear 
programming  methods  is  that  it  is  impossible  to  represent  this 
allocator  in  a  useful  format.  The  control  allocation  algorithms 
have  ddes  as  an  input  and  S  as  an  output.  In  the  pseudo 


inverse  development,  a  direct  relationship  exists  between 
inputs  and  outputs.  When  c  =  0  and  W  =  I,  Eq.  22  becomes 

6  =  Bt(BBt)-1  [ddes]  =  B^ddes  (39) 

In  this  case,  the  control  allocator  is  simply  a  gain  matrix 
given  by  B: .  This  is  useful  for  stability  or  robustness  analysis 
because  a  model  of  the  control  allocator,  namely  B:.  has  been 
determined.  When  linear  programming  control  allocators  are 
utilized,  no  such  model  exists  and  the  relationship  between 
inputs  (ddes)  and  outputs  (S)  is  more  complicated.  In  this 
case,  determining  stability  or  robustness  of  closed  loop  sys¬ 
tems  is  not  feasible.  Fortunately,  there  are  situations  where  a 
linear  programming  control  allocator  can  be  modelled.  If  no 
control  effectors  exceed  their  rate  and  position  limits,  then 
driving  the  preference  vector  in  Eq.  38  to  the  2-norm  (or 
pseudo  inverse  solution),  given  by  Eq.  39,  provides  a  method 
to  model  the  allocator.  This  is  because  the  solution  to  the 
control  allocation  problem  will  be  the  2-norm  solution  and 
hence,  the  linear  programming  control  allocator  block  can 
be  replaced  by  B:  for  use  in  stability  or  robustness  analysis. 
When  control  effectors  exceed  rate  or  position  limits,  this 
model  is  no  longer  valid.  Note  that  this  technique  can  also 
be  used  in  the  mixed  optimization  problem  described  in  the 
next  section. 

H.  Mixed  Optimization 

The  mixed  optimization  problem  is  a  combination  of 
the  control  and  error  minimization  problems.  The  mixed 
optimization  problem  is  posed  as  follows:  given  a  control 
effectiveness  matrix  B  and  a  preferred  control  vector  dp,  find 
a  vector  S  such  that 

J=  ||B<5-ddes||-H|<5-<5p||  (40) 

is  minimized,  subject  to  <5  <  S  <  S  where  v  £  ]Rn  is  a 
factor  which  weights  the  relative  importance  of  the  error  and 
control  minimization  problems.  Bodson  [7]  has  converted  the 
mixed  optimization  problem  into  a  standard  linear  program. 

I.  Affine  Control  Allocation 

A  linear  control  allocation  problem  (as  described  in  Eq.  1), 
suffers  from  the  fact  that  it  assumes  that  the  individual  entries 
in  the  control  effectiveness  (slopes  of  moment-deflection 
data)  pass  through  the  origin.  In  a  local  sense,  this  is  typically 
not  the  case  and  it  is  directly  attributable  to  nonlinear  effects 
in  the  moment-deflection  relationship.  A  more  accurate  so¬ 
lution  to  the  control  allocation  problem  is  obtained  using  an 
affine  control  allocation  problem  formulation  [5],  that  is,  one 
of  the  form:  find  S  such  that 

B<5  +  e  =  ddes  (41) 

is  minimized,  subject  to 

6  <6  <6  (42) 

where  e  is  an  intercept  term  that  provides  a  more  robust 
control  allocation  algorithm  when  the  moment-deflection 
curves  are  not  truly  linear. 
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Fig.  3.  Typical  Moment  Curve  and  Slopes  at  Operating  Conditions. 


To  begin  the  analysis  of  the  intercept  correction  term, 
consider  a  representative  one-dimensional  moment  deflection 
curve  shown  in  Fig.  3.  The  horizontal  axis  is  the  control 
effector  position,  S,  while  the  vertical  axis  is  the  moment 
produced  by  the  control  effector,  gg .  At  the  operating  points, 
(5i  and  <5-2,  the  control  effectiveness  is  given  by  the  slopes  b\ 
and  62.  while  the  intercepts  of  the  moment  axis  are  ei  and 
62-  In  the  linear  control  allocation  problem,  the  moment  axis 
intercepts  are  zero,  as  shown  by  the  slope  bi  and  62  control 
effectiveness  lines  that  pass  through  the  origin.  When  the 
control  is  operating  at  the  position  given  by  5k  (see  Fig.  4), 
the  effectiveness  is  given  by  the  tangent  to  the  moment- 
deflection  curve  at  that  operating  condition.  In  a  purely  linear 
control  allocation  problem,  this  slope  is  translated  along  the 
vertical  axis  until  it  intersects  the  origin  and,  as  shown  in 
Fig.  4,  can  have  adverse  effects  on  the  computation  of  5 
required  to  provide  the  desired  moment.  Here,  we  begin 
by  assuming  that  the  control  is  currently  operating  at  5k 
and  draw  the  tangent  to  the  moment-deflection  curve  at  that 
point  and  also  translate  this  tangent  line  to  the  origin.  When 
the  next  moment  command  is  provided,  assume  that  the 
command  is  ggdes .  If  the  control  allocation  algorithm  had 
perfect  knowledge  of  the  moment-deflection  curve,  then  the 
algorithm  could  accurately  compute  5perfect  as  the  required 
control  position  to  produce  the  desired  moment.  However, 
in  a  linear  (or  affine)  allocation  problem,  only  moment- 
deflection  slopes  are  available  to  the  allocation  algorithm.  If 
the  linear  case  is  utilized,  then  a  large  error  exists  between 
Sperfect  and  the  position  computed  using  linear  information, 
5 linear-  On  the  other  hand,  if  an  affine  representation  of 
the  moment-deflection  curve  is  utilized,  the  error  between 
Sperfect  and  Saffine  is  much  less  than  the  error  between 
Sperfect  and  5 nnear.  Hence,  this  method  can  produce  signif¬ 
icantly  more  accurate  results  for  the  computation  of  5. 

It  should  be  noted  that  the  affine  method  is  well  suited 
to  algorithms  which  obey  the  rate  limits  of  the  effectors 
and  are  implemented  in  digital  computers.  The  rate  limits 
and  digital  implementation  essentially  limit  the  distance 
a  control  effector  can  travel  in  one  timestep.  Therefore, 
even  if  operating  in  a  nonlinear  region  of  the  moment- 
deflection  curve,  the  rate  limit  restrictions  allow  the  control 


Fig.  4.  Linear  Versus  Affine  Control  Effectiveness. 


effectiveness  (tangent  of  moment-deflection  curve)  to  be  an 
accurate  representation  of  the  nonlinear  curve. 

IV.  Conclusions 

This  work  provided  an  overview  of  some  of  the  techniques 
used  to  address  the  linear  control  allocation  problem.  Meth¬ 
ods  to  reduce  the  dimension  of  the  control  effector  space, 
such  as,  explicit  ganging  and  daisy  chaining,  were  presented. 
More  sophisticated  algorithms,  that  take  into  account  control 
effector  rate  and  position  limits,  were  also  discussed.  In  these 
cases,  linear  programming  has  proved  to  be  a  viable  method 
for  solving  online  control  allocation  problems. 
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